Possibility of long-range order in clean mesoscopic cylinders 



M. Lisowski, E. Zipper, and M. Stebelski 
Institute of Physics, University of Silesia, ul. Uniwersytecka 4, 40-007 Katowice, Poland 

A microscopic Hamiltonian of the magnetostatic interaction is discussed. This long-range inter- 
action can play an important role in mesoscopic systems leading to an ordered ground state. 

The self-consistent mean field approximation of the magnetostatic interaction is performed to give 
an effective Hamiltonian from which the spontaneous, self-sustaining currents can be obtained. 

To go beyond the mean field approximation the mean square fluctuation of the total momentum 
is calculated and its influence on self-sustaining currents in mesoscopic cylinders with quasi-lD 

Qand quasi-2D conduction is considered. Then, by the use of the microscopic Hamiltonian of the 
magnetostatic interaction for a set of stacked rings, the problem of long-range order is discussed. 
OO ' The temperature T* below which the system is in an ordered state is determined. 
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I. INTRODUCTION 



One of the most exciting areas of physics is the study of mesoscopic electronic systems, i.e. metal or semiconductor 
samples which are sufficiently small and at sufficiently low temperature, such that inelastic elecijijon-phonon scattering 
is reduced and the electron propagates as a phase coherent wave throughout the entire sampleEl. 

Recently in a series of papersHa we discussed a possibility of spontaneous persistent currents in relatively clean 
(ballistic regime) metallic or semiconducting systems of cylindrical geometry. 
'"^ , It was shown in the mean field approximation (MFA) that the inclusion of the magnetostatic interaction among 
' electrons can lead to an ordered ground state with spontaneous self-sustaining orbital currents which run without 
O ' support of external magnetic field. 

I ^ I In our investigations we considered a collection of many mesoscopic rings with a thickness d R, {R is the radius 

of the ring) stacked along z axis and the three-dimensional (3D) mesoscopic cylinder of very small thicknessQ. 
T-H , In this paper we want to give some justification to the calculation mentioned above, because there was no microscopic 
^ ' theory of this phenomenon till now. We will examine a microscopic Hamiltonian for the magnetostatic (current- 
current) interaction and show that the self-consistent MFA of it gives the effective Hamiltonian H^^^ leading to 
self-sustaining currents. We also show that magnetostatic interaction is long-ranged and therefore the criteria of the 
MFA are met. 

The orbital magnetic interaction and its static version, the magnetostatic coupling, has been discussed previously by 
00 ' Pines and Noziereaj and has been shown to be small in macroscopic metallic samples. However this interaction should 
On I be reexamined in mesoscopic systems where the presence of energy gaps in the energy spectrum changes qualitatively 
its physical properties leading e.g. to persistent currents driven by the static magnetic flux (j) at low temperaturesQ. 
j3 . Persistent currents create orbital magnetic moments and their interaction can lead tO|jttteresting coherent collective 
phenomena which bear some resemblance to ferromagnetism and to superconductivityax 

To go beyond the MFA and discuss fluctuations we follow the ideas developed by Blochcl. He discussed the problem 
of quantum coherence in a macroscopic, metallic ring manifested e.g. by the thermodynamically stable flux trapping 
Q at zero external magnetic field. He covers in his paper different long-range characteristics of the normal and the 
Q , superconductive state of a metal. 

^ ' Irrespective of the specific dynamical properties of the systems he shows that the general criteria for fiux trapping 
• i-H , are closely related to the mean square fluctuation of the total momentum and depend strongly on the dimensionality 
' of the system. 

^ We will show, using Bloch's formulaes that if we reduce the dimensions of a macroscopic cylinder made of a normal 

■ - - ■ metal or semiconductor to mesoscopic dimensions, the system exhibits coherent properties absent in macroscopic 
samples. We will formulate the criteria under which flux trapping can be obtained in mesoscopic cylinders with 
quasi-lD and quasi-2D conduction. In particular we will discuss the mean square fluctuation of the total momentum - 
the decisive quantity for the characterization of the properties of the system. We show that it is smaller in mesoscopic 
systems than in the corresponding macroscopic ones thus favouring quantum coherence. 

Finally, using the microscopic Hamiltonian we discuss the possibility of the long-range order in a mesoscopic cylinder 
made of a set of mesoscopic rings. We calculate the correlation length and the characteristic temperature under which 
the system is in a magnetically ordered state. 
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II. MAGNETOSTATIC INTERACTION 



The general formula for the magnetostatic interaction is of the form: 



Hmgt=-— drdr- — , (1) 



in 



where J.{r) = ep{r)/me, J.{r) is the current density, p{r) is the momentum of an electron. 

Let us assume that the currents run in a set of mesoscopic rings of small thickness deposited along z axis. We 
can write 

= (2) 

m/l •^'^^ 



where Cm is given by a parametric equation for an electron going around the circumference of a ring, ^ — ^ (s), s 
is the coordinate along the circumference of the ring, Im is the current in the m-th ring. 
We obtain 



M. 



Hmgt 2 

m/l m'/l 

where 



mm — . 

An 



c„ JC, 



—7n —m' 



^mm' — ^m'm- (4) 



Thus we have obtained the interaction Hamiltonian of the currents from different rings. In this derivation we 
neglected the selfinductance effects in each single ring. It can easily be seen that they are small. 

The interaction constant C„im' depends on the sample geometry; here it has to be calculated for the rings deposited 
along z axis at|jdistance Zmm' = Zm — ^m' • 



The result isO: 



^mm' — /^O-^ 

where 



K - -^E 



(5) 



- - (6) 



mm' 



^ ^mm' 



Jo fl -/,2cin2.^l/2 J, V ; 



(7) 



K and E are the elliptical integrals of the / and // kind respectively. 



The z dependence of the coupling constant £ is presented in Fig. 1. We see that the interaction (||) is a long- 
range interaction. For small z it falls down slowly proportionally to fioRln (R/z — 2), for large z it falls down faster 
proportionally to 1/z'^. The interaction constant depends only on R and on the relative distance of the centres of the 
rings. 

As the currents Im can run only in the clockwise or anticlockwise direction the Hamiltonian (H) has the form of the 
Ising Hamiltonian. It can be also expressed via the momenta pm, Pm' from different rings. 
For the ring geometry we get 
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Pm — ^ ^ Pnm ^ 

n/1 

where N^' is a number of conducting electrons in a single ring. 
Hmgt given by Eq. (0) can be rewritten in the form: 



o2 



Hmgt — 2m^ ^ y ^ ^ gmm'PmPm' , (9) 
m/1 m'/l 

If we add to Hmgt the kinetic energy term and if we assume that the external magnetic field parallel to the z axis 
is applied to the system, we obtain Hamiltonian H: 

^ = ^Y.Y.Pnm- l^Y^Y. 9mm-PmPm' , (10) 
m/ln/1 ^m/lm'/l 

where Pnm — Pnm ~ 6^e; Pnm ^hc momentum of the n-th electron in m-th ring. is the vector potential of an 
external magnetic field which points in the x direction measured along the ring. 

The first term in Eq. (^0|) represents the kinetic energy of electrons in the external magnetic field, the second term 
represents the orbital magnetic interaction; due to the negative sign it favours states with large total momentum in 
competition with the kinetic energy. 

III. MEAN FIELD APPROXIMATION OF MAGNETOSTATIC COUPLING 

We have seen in chapter 2 that the magnetostatic coupling is a long-range interaction. This indicates in particular 
that thermodynamic fluctuations of the current will be strongly supressedt2l. 

Let's perform a self-consistent MFA of the interaction (|^), such approximation is known to be good for a long-range 
interaction. 

Hmgt = -l^^Pm^gmm'Pm' = ' ^^'^P^nAl{Zm) , (11) 

^ m m' m 

e 



^li^m) — T: ^ ] gmm'Pm' , (12) 



2m 

m' 

BC^B (C) + {B)C- (C) {B) 

Hmgt - -TT^ [2pm {Al{Zm)) ^ (Pm) {Al{Zm))] , (13) 

m 

where the first term in Eq. ( p^ has been obtained by use of the symmetry relation (^, {Aj{zm)) ~ 
(e/2me)^ , gmm' {Pm') = Ai where we assumed that (pm) = (p) is the same for all rings, and hence from Eq. 
(i (Im) ^ (/). 

Assuming that our stack of rings forms a long cylinder of length I we can calculate the vector potential Aj, We 
get 

AJ^^,oR^^. (14) 
Calculating the current (/) with a total vector potential A = Ae + Aj we get 
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Inserting (|5|) into (jlj) we obtain the self-consistent equation for Aj 



Aj = Hp) - N^eAj) , (16) 



from which we get 

^ 1 + riN^M, 



= T-^hwrrMz (p) , (17) 



where 



Inserting (p) from Eq. ( p^ into Eq. dl^ ) we obtain the Hamiltonian H (Eq.(|l0|)) in the self-consistent mean field 
approximation: 

'^-•"-iEE^»™-»V + |, (18) 

m/1 n/1 

where C = fioTrR^{VP - - R)/l^, (pi = 2'kRAi. „ 

The Hamiltonian (|l^) was the basis of our previous investigations^ of spontaneous self-sustaining currents (see Eq. 
(^), (|^))- Its derivation from the long-range current-current interaction serves as a justification of the use of the 
Hamiltonian ( [l8| ) to investigate magnetic properties of mesoscopic systems. 

It has-|been generally believed that the MFA should work well for long-range forces. However, it has been shown 
in Ref.Eil that the above statement is correct if an additional condition is fulfilled. The authors defined there the 
quantity 5: 



\2 



^y~T^ , (19) 

mm' 



where Lmm' is the interaction constant. They proved that MFA is correct if ^ 1. 

We have calculated S with Lmm' given by Eq. (||) for the following set of parameters: h = Zm,m+i = lOA, 
R = 5000A, ~ 10^. We have obtained S ~ 10^; it means that the MFA should work well in the case considered 
by us. 



IV. SUMMARY OF THE BLOCK'S RESULTS 



We ar&|going now to investigate the possibility of flux trapping in mesoscopic cylinders using the ideas developed 
by BlochB. The general criteria for flux trapping are closely related to mean square fluctuations and give a natural 
way to describe the system in terms of the two-fluid model. At first we briefly recall Bloch's results. 

Let us consider a system of N particles with mass m and charge e, contained in a ring with radius R and radial 
width d <^ R. We assume that the magnetic field parallel to z axis is caused by a current around the ring, so the 
vector potential A = 0/(27r_R) points in the x direction. 

The total momentum in the x direction is given by 

N 

P.=T.Pn' (20) 

n/1 

where p° represents the momentum of n-th particle {n — 1,2, N). 
The Hamiltonian of the particles is of the form: 

H^^-^^-J^^H', (21) 
2Nme 
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where H' contains the kinetic energy of the motion in the y and z direction and of the relative motion in the x direction 
as weh as any additional terms which arise from interactions and characterize the specific dynamical properties of the 
system. P 

From the symmetries and periodic boundary conditionsQ we get the eigenvalues of Px'- 

P^n^ = {Nu + fi)^, (22) 

n, v are arbitrary integers, fi is likewise an integer such that 

N N 

< u < — . 

2 ^ - 2 

Eigenenergies of the system are given by 

' "2^^ + 

q represents the system of additional quantum numbers necessary in addition to P {v, fi) in order to fully characterize 
the state of the system; 0' — (t)/(j)o, <j>o — h/e. 



Using Eq. (23) we can calculate the free energy from the particles. The flux dependent part of the total free energy 

is 

F{cp') = F,{cp')+F2{cp'), (24) 
where Fi{(f)') is a periodic function of with period 1, 

Fi{<j>') = -kBT\nZi{<j>'), (25) 



■^i(<?^') = \It77 I l + 2^age"^cos27r50' ] , (26) 

9/1 
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^z^e-2"fM/JV^ (27) 



7 ~ I (2meR^kBTy, is the statistical weight of the state E'^^ , 

= 1. > 0- (28) 



= (29) 

F2{4>') is the energy stored in the magnetic field, C is the self-inductance of the ring. 

Thermodynamically stable flux trapping is determined by those values of (j)' for which F{(j)') has a minimum. To 
achieve it a strong variation of Fi{(j)') is necessary to prevent the dominance of -F2 (</>') which has minimum at 0' = 0. 

Let us consider now three special cases which elucidate how different states of matter can be described in this 
model. The information about specific properties of the system are contained in the quantities z^. 

A. z^ is independent of /i , i.e. z^ — 1/N. 

After some algebra we arrive at 

, / oc „ 

e — cos27r.giV(/)' . (30) 

The general periodicity in 0' = 1 is accompanied here by a far shorter period (j)' = 1/N what is vanishingly small for 
large N. Besides the amplitude of the oscillation is small at any realistic temperature leading to a very small variation 
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of Fi{(j)') and hence to the absence of stable flux trapping. This situation corresponds to the case where the system 
exhibits no long-range order and is characteristic of the normal state of a metal. 



B. = Sf^Q 



From Eq. (|26|) wc find: 




(t3) 



Z^{(j)')= —11 + 2} e~T^cos27rg0' , (31) 



which has the periodicity with = 1. However in order to get a stable flux trapping for 0' 7^ we need 



what is satisfied at low temperatures T. Then fi (</>') can well dominate the part F2{4>') in Eq. (p4[). The pronounced 
minima of F{(j)') occur &t 4>' = v and are equivalent with stable flux trapping. This case corresponds to a condensed 
ideal Bose gas where all particles have momentum p = vh/R. 

The partition function is of the form: 

I I °° 2 \ 

^1(0') = ^^ ( l + 2^e"^cos47rg<^'j . (33) 

This expression has the periodicity in 0' with (f)' = ^ i.e. with = h/2e. This case exhibits the property that N/2 pairs 
have all the same momentum p — vh/R and corresponds to the long-range order characteristic of a superconductor 
at temperatures in which the condition (52) is satisfied. 



To consider other, more general cases one assumes that the total momentum P can admit the values given by Eq. 
( p^ ) with /i 7^ 0, but with a sharp maximum around fi = 0. The formula for Zi{(j)') has then the form: 

I 1^2 \ 

Zi((/)') = ^^ ( l + 2^e^^^cos27rg0'j , (34) 



where 



P= \j , 0<p<l, (35) 



{APfj is the mean square fluctuation of momentum P around the set of values (l)'Nh/R, where (/)' is an integer 
multiple of the period, p can be called a relative fluctuation. 

{(APf) . ^ ((A,)^) . -mefeBT^P.^, (36) 

a 

where (^(AF)^^ = (P^) - {Pf ; /(p„) is the mean number of particles with momentum pa = h {a — (j)') / R, a = 
0,±1,±2,... 

Eq. (^) covers the special cases discussed before. The case A corresponds to the situation where all values of ^ 
are equally probable what leads to p ^ 1; Zi{(j)') becomes then independent of </>' with the exclusion of flux trapping. 
Cases B and C correspond to p = and thus describe the systems in the coherent state with no fluctuations. 

The intermediate case: 

< ((AF)^J> < Nm^kBT (37) 

or equivalently < p < 1 still results in a pronounced variation of Fi[(j)') as long as 

A^7(l - p) > 1. (38) 



It means that stable flux trapping can be expected as soon as ( (AP) ^ is found to be smaller than the maximal value 
of mean square fluctuation of P. 
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V. FLUX TRAPPED IN MESOSCOPIC CYLINDERS 



We will perform now, using Bloch's formalism, some model calculations of flux trapped in mesoscopic hollow 
cylinders of radius R, length I and wall thickness d (d R) made of a normal metal or semiconductor. Such cylinders 
can be treated as a multichannel system with channels in the length and Mr channels in the thickness of the 
cylinder (M = M^Mr), with the total number of conducting electrons N — N^M. We will study the systems ■^ptbi 
quasi-lD and quasi-2D conduction. It is known that coherent response in mesoscopic cylinders can be obtainedlla'cl 
for systems with large phase correlation of currents from different channels which is related to the shape of the Fermi 
Surface (FS). The most favourable situation is for systems with quasi- ID conduction, i.e. with a flat FS parallel 
to axes of the wave vector k. There exists then a perfect corelaticm among the channel currents and the magnetic 
response is the strongest. Such FS can be obtained in bcc crystalsllj, in low dimensional organic conductors with the 
overlap of the orbitals mainly in one direction, and when a cylinder is made of a set of M quasi-lD rings stacked 
along z axis by e.g. lithographic method. Such multiple quantum chains can be mapped into a system with the flat 
(rectangular) FS. 

In the case of quasi-2D conduction our cylinder can be constructed from a set of M,. 2D coaxial cylinders (e.g. 
multiwall carbon nanotubes and cylinders nwde of a material with layered structure) . 
To simulate different shapes of the 2D FSailil we can use the equation: 

kl = k'P^ + k^p^, (39) 

where u is an integer number. 

For u = 2 we get the circular FS, for u > 12 the rectangular one and for 2 < u < 12 the rectangular FS 
with rounded corners. For the circular FS the currents from different channels add almost without correlation, the 
correlation increases with increasing the curvature of the FS, i.e. with increasing u. 

Let us consider at first a set of M quasi-lD rings stacked along z axis (or in general a cylinder with quasi-lD 
conduction). Let us assume that the magnetic field is caused by persistent currents from all rings. Thus we meet 
conditions from the Bloch's paper and our mean field Hamiltonian H^'^^ (Eq. (|l|)) leads to the free energy given by 
Eq. (H) with 

F^{ci,')^MF^{<ly'), (40) 

where F^{(f)') is the free energy of a single ring. 

In general an external magnetic flux (p^ parallel to the z axis can be also applied to the system but we are mainly 
interested in the self-sustaining flux at 0e = 0. 

The necessary condition for coherent behaviour has the form: 

iV^^^»l, (41) 

where Aq — / {2meR^)^ Aq is the quantum size energy gap at the FS. 

Let us calculate the momentum P defined as in Eq. ( p^ ) for three different model cases. 
1. If the number of electrons in each ring is odd we find: 

p^Nv— ior [v--] <(j)' <{u+-], i/^0,1,... (42) 



R V 27 " V 2, 

what corresponds to the case B with z^ — <5^o p = 0. 

Calculating the minimum of Ihe total free energy (Eq. (|24|)) dX (pe — Q and at T ^ Ao/ks we get the value of the 
flux trapped (/)'* in the cylinderB: 



^odd 



1 + K 



(43) 



where 



— 2 — and M<-(l + K,-^] 
c2r AT I I — o V ' ) 



e^CN ' ' - 2 

Eq. ( p3| ) reveals the influence of the finite size effect - the fiux trapped is quantized in units less than 000. For 
macroscopic samples k ^ and (/>q^^ = ly. 
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2. If the number of electrons in each ring is even, the momentum P is 













H) 









<0'<(j^ + J), 1^ = 0,1,... (44) 



2 



It corresponds to Zn = S n_ and stable values of trapped flux are 



even 



1 



< 1(1 + (45) 



For macroscopic samples 4''lvm — ^ ~^ \ what corresponds to stable minima of the free energy at half-integral values 
of 00- 

We see that a system under consideration exhibits thermodynamically stable persistent currents and flux trapped 
at temperatures T <C Aq/Zcb. This condition is easily satisfied for mesoscopic rings at T < \K (e.Sr-|for R ^ \ /im, 
^a/ks ^ 12.5i4r), it is however unrealistic for macroscopic samples (for R ~1 cm, Ao/fcs ~ 1Q~^K]^. 

Our treatment is not only for identical rings. We have performed the calculations for a set of rings in which the 
number of electrons changes in the range iV« = ± AiV-^, AiV-^ ==10. We have considered two kinds of 
changes: 

i) N^^ fluctuates from iV^ to iV^ ± 2n, n — 1, 2, 3, 4, 5 - the influence of such fluctuation on the persistent current 
is very small. 

a) changes from to ± (2n — 1) - the influence of this kind of changes on the current is pretty large (see 
below). 

In the presented paper we study mainly the systems with the diamagnetic reaction on small magnetic flux (jj. The 
rings with odd give a diamagnetic current whereas those with even give a paramagnetic current. We have 
found that the diamagnetic reaction and trapped flux can still be obtained if about 20% of rings carry an even number 
of conducting electrons. 

3. Finally, we may also consider the model case where roughly a half of the rings has an even number of electrons 
and a half of the rings has an odd number of electrons. It corresponds to the case C where = ^ (^^^^o + ^f^Jij ^nd 
the minima of the free energy occur with the twice smaller period. r. 

We have not discussed this situation here in details because the total current in this case is paramagneticB. 

We are going now beyond the MFA and consider fluctuations around the set of P values discussed above. We 
assume that the total momentum P in a cylinder given by Eq. ( |2^ ) can admit the values with /i ^ and/or /i ^ N/2 
but with a maximum at = and/or fi = N/2. This assumption leads to the partition function Zi((j)') given by 
Eq. (^^ and the criteria for quantum coherence manifesting themselves in the flux trapping are ultimately related to 
the magnitude of the mean square fluctuation of momentum P. The condition for coherent behaviour takes then the 
form: 

Ao(l-p) , . 

» 1> (46) 

and we see that the fluctuations decrease the energy gap. Eq. ( |3^ ) covers also as a special case a situation characteristic 
of a normal state of a metal (no flux trapping) in which all values of fi can be found, with equal probability, in the 
formula for P. 

Stable flux trapping can be expected as soon as (^(AP)^^ is found to be a small fraction below NmeksT, called 

by Bloch the equipartition value. In the following we will calculate {^(AP)^^ from Eq. (|3^ ) and relate its magnitude 

to the coherent behaviour of the sample. 

The mean number of electrons with momentum Parm fiPam) is given by the Fermi-Dirac distribution function: 

/(p„„0 = /(-B„,„(0)) = (47) 



where /ic is calculated from the condition: 



N = MrJ2 E /(^-m(</>)). (48) 

m/l a/0,±l 
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The electron energy eigenvalues, calculated by the use of periodic boundary conditions in the x direction and cyclic 
boundary conditions in the z direction: 



2m„ 



— a — eA] +h k^{m) 



is the total flux contained in the cylinder, (jjj — CI{4)), where 

Ah ±00 



(50) 



^^^^ = 9^^^^'- E E - <i^')f{E^mm, (51) 

m/l a/0,±l 



and after expansion in the Fourier series! 



m = ^^Mr E 5] fcp^ (m) ^ /^fr'ix cos(2^5i?fcF. (m)) sin(2^#'), 



m/l 5/1 



l-exp(-^) 



where a is the lattice constant, and according to Eq. ( ^ ) kp^{m) = hp [1 ~ {kz{Tn)/kF)"]^^" , kz{m) — rmr/l, 
m= 1,2,...,A4 



Eqs. ( pOD and (51) form a set of self-consistent equations for the current. The question of existence of self-sustaining, 
persistent currents is reduced to the problem whether these equations have stable, nonvanishing solutions at (f>e — 0. 

Notice that the dispertion relation (^) is modified by the presence of the flux (pi coming from the currents and has 
to be calculated in a self-consistent way. We show below that 0/ produces a dynamic gap in the system and therefore 
increases coherence. 

The values of (|(AP)^^ for different shapes of the FS and for R ~ 10'*A/(27r), = 10000, Mr = 100 at T = 15K 
are presented in Table 1. 

For the rectangular FS, corresponding to quasi-lD conduction, we assumed that it lies in the middle of an energy 
gap Aq for an electron going along the circumference of the cylinder. We see that the magnitude of (^(AP)^^ (or the 
corresponding relative fluctuation p) decreases with increasing the curvature of the FS. The difference between the 
normal and the coherent state of a mesoscopic cylinder is reflected in the magnitude of ^(AP)^ 

In the ideal, limiting case corresponding to (^(AP)^^ — (case 1 on page 0) the system is fully coherent i.e. 
N — Nc, where Nc is the number of electrons in a coherent state. 

Finite values of ^(AP)^^ can be interpretated by use of the two-fluid modeflS as being proportional to Nn = N—Nc, 

Nn is the number of particles in the normal state. We can write {^(AP)^^ — NniriekBT. 

The maximal values of (|(AP)^^ given by Eq. (|3^) will be obtained for macroscopic values of R where the 

replacement of the sums over a by integrals is permitted. This gives us (^(AP)^^ — NmeksT or N = Nn and the 

system is in a normal phase. However for mesoscopic cylinders such replacement is not allowed and the presence of 
finite size energy gaps leads to Nn < N what means that a part of the electrons is in a coherent state. The presence 
of coherent electrons results in persistent currents, the amplitude of which depends on the shape of the FS. 

Persistent currents driven by an external flux 0e vanish if we switch the external field off. However the presence of 
the flux coming from the magnetostatic interaction can lead to persistent self-sustaining currents or in other words 
to flux trapped. In general persistent currents can be paramagnetic or diamagnetic. Paramagnetic self-sustainiiig 
currents correspond to spontaneous current^, and diamagnetic self-sustaining currents correspond to flux trapping. 
In the following we will study mainly the diamagnetic solutions. 

In order to discuss the influence of on the coherent properties of mesoscopic cylinders we calculate an energy 
gap at the FS for electrons going around the circumference of the cylinder (and for (p < 0o/2). Using Eq. ( ^ ) we 
find: 

Ap = E^^ + i^ra - E^^^ra = aJ I - 20^ + 2^) . (52) 
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Ap contains a term 



£ I/I 

A, ^ Ao^, (53) 

Ad is the dynamic part of an energy gap which should increase coherence in the sample. 

That this is really the case we can see from the comparison of the first and second column in Table 1. We see that 

the mean square fluctuation ^(AP)'^^ when calculated with the dispersion relation pertinent for normal electrons 



namely with £am = [{ha / R — eAe)'^ + h kf {m)]/{2me) is larger than ^(AP) j calculated with the dispertion relation 

(^), modified by the presence of self-consistent flux Thus the presence of the magnetostatic coupling decreases 
fluctuations and the number of normal electrons in the sample. This is also seen in Fig. 2 where the temperature 
smearing of the distribution function calculated with the self-consistent flux is smaller than that calculated with 0e- 

We also checked that, as should be expected, the mean square fluctuation {^(AP)^^ decreases with increasing the 

number of interacting channels. _ 

This way of analysis bears some resemblance to the two-fluid description of a superconductonlZI. For free electrons 
with the dispersion relation the electron density can be written as 

k — 

Defining similarly the density of normal electrons in a superconductor, with a dispersion relation Ek = V^k + 

(55) 

and the density of super-electrons as 

nc = n- Un, (56) 

we get the two-fluid description. 

We are going to discuss now the influence of fluctuations on self-sustaining currents. In Fig. 3 we present the 
currents given by Eqs. (|50|), ( pT] ) at 0e = 0, for different shapes of the FS and for different values of p. The 
fluctuations decrease the current, but self-consistent solutions can still be obtained for Fermi surfaces with flat regions 
{u — 6 and u — 12). However, as should be expected, the value of self-consistent current (flux) is smaller for the case 
with fluctuations included. For u = 2 we do not get flux trapping because the number of coherent electrons is too 
small. 



2h^ 



k 



Finally, to get more insight into the properties of our system we discuss, using the microscopic Hamiltonian (|^), 
the possibility of long-range order in a cylinder made of a set of mesoscopic rings. 

In the literature one finds the statements that phase transitions and long-range order are impossible in quasi- ID 
systems. It is true for systems in the thermodynamic limit with short-range interactions. However, one can look for 
the conditions for an ordered state for large but finite number of interacting entiticst^. 

Let us consider a set of rings described by the Ising-like Hamiltonian (^). In the ground state all currents run 
parallel. Let us construct a new configuration by reversing the direction of L currents {L <C Mz) in s (1 < s <C Mz) 
different places in the chain of rings. The energy change will be denoted by sA£'i(Mz), and the change in the free 
energy is 

AP = sAEl{Mz) - TAS = s \AEl{,Mz) - ksTlnMz] . (57) 
If AP > then the ordered configuration is stable. This condition is equivalent to 

= e^^-/'^^^ > Mz, (58) 

where has the sens of a (dimensionless) correlation range. 

If Eq. ( |58| ) is fulfilled the system is ordered in the sense that the correlations extend over all its length. 

Let us discuss the possibility of long-range order for the case considered by us. The energy AEl{Mz) is of the form: 
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^ / „/l m-m'/L + l 

The L dependence of AEl{Mz) is presented in Fig. 4. We see that AEl{Mz) increases with L and decreases with 
b = Zm,m+i = Zm+i — Zm- The smahcst energy change is obtained when reversing a direction of a single current in 
several different places. In Fig. 5 we present the dependence of Ai?i. We see that Ai?i increases with for 
small Mz and then saturates. 

We are in position now to calculate the temperature T* at which the crossover from an ordered to disordered state 
occurs: 

The calculations performed for the following set of parameters: — 10'*, — 1, R — 5000A, b = lOA gave us 
T* ^ O.IAK. It means that at T < T* the system exhibits a long-range order in the sense that the correlation 
range is longer than the sample size. The temperature Tc calculated in the MFA for the above set of parameters is 
Tc ~ 0.216KB 

We see that T* obtained by the use of the Ising model is of the same order as Tc obtained with the MFA. It means 
that long-range interactions encountered in our system strongly supress fluctuations. If we assume that each ring has 
a small number of transverse channels Mr then e.g. for Mr ^ 3, T* ~ IK. 

It has to be stress e.g. that finite value of T* for our system is a finite size effect. Indeed, as we can see from Eq. 
® and Fig. 5 T* ^ for M^^oo. 



VI. DISCUSSION AND CONCLUSIONS 



In the presented paper we discussed the magnetostatic coupling of electrons in mesoscopic systems and its approx- 
imations. This interaction is known to be weak in macroscopic samples, however it seems it can play an important 
role in mesoscopic samples due to very peculiar properties of mesoscopic systems in the magnetic field. 

It is possible to induce persistent currents (or in other words orbital magnetic moments) in mesoscopic ring by the 
static magnetic field. The magnetic interaction of orbital magnetic-jnoments can lead to magnetically ordered ground 
state. This possibility has been discussed in a number of papersfl'Ell using mean field approximation which states that 
each electron moves in an external magnetic field and the field coming from all currents in a system. The obtained 
two self-consistent equations for the current can lead to spontaneous self-sustaining current at zero external field. We 
neglected here the Zeeman energy of the electron spins becaarse it turns out to be very small compared to the orbital 
energies. The infiuence of spins has been discussed in Ref.E£l. 

In this paper we have presented the microscopic Hamiltonian which is responsible for the internal magnetic field 
- the magnetostatic (current-current) interaction. The strength of the interaction depends strongly on the sample 
geometry. For the stack of mesoscopic rings deposited along certain axis we get the long-range interaction with the 
coupling constant depending only on the radii of the rings and on the relative distance of its centers. 

We have shown that the self-consistent MFA of the current-current interaction, gives the effective Hamiltonian 
jjMF leading to self-sustaining currents. Its derivation from the long-range interaction serves as a justification of 
the use of the Hamiltonian (^8|) to investigate magnetic properties of mesoscopic systems. The MFA is known to be 
the best for systems with long-range forces (if the condition (^9|) is fulfilled), thus we should expect that it leads to 
reasonable results in the considered case. 

To obtain the full Hamiltonian describing our system one should add to Hamiltonian given by Eq. (|l^) the Coulomb 
interaction and the interaction with impurity jpntential. It was recently shown that the Coulomb inteeaction does not 
infiuence persistent currents in clean systemsc2l, whereas it enhances the current in diffusive regimecil. In a work by 
M. Pascaud and G. MontambauxEj the experiments which permit to test the role of Coulomb interaction have been 
suggested. 

In the model calculations presented in this paper we did not consider the effect of impurities in order not to obscure 
the whole subject with too many details. The influence of disorder on self-sustaining currents has been analysed 
in Ref.EI. We found that disorder decreases persistent currents but self-sustaining currents can still be obtained for 
relatively clean samples (ballistic regime). 

To go beyond the MFA we have considered the influence of fluctuations, calculated by the use of Eq. (|3^), on the 
properties of a mesoscopic cylinder. We have shown that these fluctuations are smaller in mesoscopic systems than 
in macroscopic ones because of the quantum size energy gaps. On the top of it the magnetostatic coupling modifies a 
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dispertion relation and creates a dynamic gap what leads to further reduction of the fluctuations. Thus in mesoscopic 
systems coherent and normal electrons coexist and the system can be described by the two-fluid model where the 
fluctuations are proportional to the amount of normal electrons. Self-sustaining currents run by coherent electrons 
survive fluctuations in systems with FS having flat regions, however their magnitude is reduced. 

Having the microscopic Hamiltonian for electrons interacting by magnetostatic coupling (Eq. (^)) we have dis- 
cussed, for a set of stacked rings, the possibility of long-range order. Phase transitions and long-range order are 
possible in the strict mathematical sense only for systems in the thermodynamic limit. For finite (but still large) 
systems discontinuities of thermodynamic quantities and infinite range correlations are not necessary. "Discontinu- 
ities" have finite widths jmd correlation ranges may be as large as the system itself, regardless of the behaviour in 
the thermodynamic limiillS. What is more finite systems can show interesting effects which will be wiped away in the 
thermodynamic limit. This is the situation in the presented paper where at temperatures T < T* ^ O.IK — IK the 
set of mesoscopic rings (or in general the mesoscopic cylinder with quasi-lD conduction) can exhibit the long-range 
order, but T* ^ for ^ cx3. 

The statistical properties of the system following from Eq. (0) will be presented in a subsequent paper. 
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FIG. 1. The interaction constant jC as a. function of distance z between the ring centres. 



FIG. 2. The Fermi-Dirac distribution function f{Ec,) versus energy Ea, for the cylinder made of a set of quasi- ID mesoscopic 
rings, in the magnetic flux (ji^ and (f) = (j)^ + ifi respectively. 0/ has been calculated in the self-consistent way for parameters 
as in the figure. 
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FIG. 3. Persistent currents / as a function of magnetic flux 0/0o in quasi-2D mesoscopic cylinders with different shapes of 
the Fermi surfaces, with and without fluctuations. Self-sustaining currents Ig. 



FIG. 4. The energy change AEl/W, where W = \ y^^^^j i ^ ^ function of the number of magnetic moments L for 
different values of 6 s Zm,m+i- 



FIG. 5. The energy change AiJi and the temperature T* as a function of the number of channels Mz in a single cylinder. 
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TABLE I. The relative mean square fluctuations of the total momentum P, p = < (AP)^ > /NmeksT , in a system made of 
2D coaxial cylinders with iV=10001xlO® conducting electrons (lattice constant a=lA, 6=5A) at temperature T=15K, in the 
magnetic flux (f>e and (p = 4>e + (pi respectively, for different shapes of the Fermi surface. 4>i has been calculated for each 0e 
in the self-consistent way, assuming the number of interacting coaxial cylinders Mr = 100. Nn represent the total number of 
"normal" electrons in the system. The number of coherent electrons is = N — N„. 



14 




10-9 10-8 10-7 10-6 10-5 10-4 10-3 10-2 

z[m] 

Fig. 1 . [M. Lisowski, E. Zipper, and M. Stebelski - On the possibility of long-range order ...] 




6.020x10-18 6.022x10-18 6.024x10-18 6.026x10-18 6.028x10-18 6.030x10 

2. [M. Lisowski, E. Zipper, and M. Stebelski - On the possibility of long-range order ...] 



0.3 



0.2- 



0.1 - 



0.0 



■0.1 - 



■0.2- 



■0.3 



1 u=12, p=0.000 

2 u=12, p=0.069 

3 u=6, p=0.000 

4 u=6, p=0.070 

5 u=2, p=0.000 

6 u=2, p=0.985 



'/// 



f/'l 



N = 10001 x106 
a = 1 A, b = 5 A 

T= 15 K 



— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.C 





3. [M. Lisowski, E. Zipper, and M. Stebelski - On the possibility of long-range order ...] 



2.5x105- 




200 400 600 800 1000 

L 

Fig. 4. [M. Lisowski, E. Zipper, and M. Stebelski - On the possibility of long-range order ...] 




10-2 I — — — — — I 

101 102 103 104 105 106 

M 

z 

Fig. 5. [M. Lisowski, E. Zipper, and M. Stebelski - On the possibility of long-range order ...] 



